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Despite its successes, the standard model (SM) of 
particle physics leaves many important questions unan- 
swered. For example the SM does not provide a candi- 
date for dark matter (DM). Direct detection experiments 
such as DAMA [1], CoGeNT 0,i3,l and CRESST 4] have 
reported signals suggestive of DM with mass in the few 
GeV/c^ range, and with coupling to the SM sector of a 
strength enabling its detection at collider experiments. 
Many beyond-the-SM theories predict DM candidates to 
include such coupling between the DM and SM sectors. 

In the framework of effective field theories, production 
of a DM particle (D) in association with a single top 
quark at hadron colliders has been recently studied [sj-LZj . 
Here, we denote the final state containing one top quark 
and dark matter as monotop. Such studies are also in- 
spired by the models of monojet produced in association 
with missing energy used to probe gravitons Q . Mono- 
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top DM production is described by a set of Lagrangians 
incorporating all possible types of DM particles (scalar, 
vector, fermion, etc.) and their flavor-violating interac- 
tions with quarks |9l4l3| . This effective theory can inclu- 
sively describe many beyond-the-SM models. For exam- 
ple, new physics processes with a monotop final state can 
also arise from the decay of a supersymmetric squark into 
a neutralino and a top quark, from the decay of a vector 
leptoquark into a massless neutrino and a top quark, or 
through flavor-changing neutral interactions with a new 
vector state escaping detection. 

In the SM, top quarks are primarily produced in pairs 
at particle colliders. They can also be produced singly via 
weak interactions, resulting in a final state consisting of 
a single top quark with additional lighter-flavor quarks. 
SM single top-quark events in the missing energy plus 
jets channel have been studied within the standard-model 
hypothesis . As the beyond-the-SM monotop theory 
predicts production of a single top quark in association 
with a DM particle, the published SM single top-quark 
results do not provide any conclusive information on the 
existence of monotop. In addition, searches for the asso- 
ciated production of top quarks with DM particles have 
only been performed in the context of events contain- 
ing a pair of top quarks Therefore, a dedicated 
search for monotops produced in colliders is needed, as 
the observation of monotops would be a clear sign of new 
physics. In this Letter, we report the first direct search 
for monotop signatures at particle colliders, assuming the 
top quark to be produced through flavor-changing inter- 
actions of up and top quarks, in association with a DM 
candidate D. We assume that the D particle has a mass 
in the range of — 150 GeV/c^; we do not consider decays 
of the D particle to up and top quarks in a higher mass 
range. 

The top quark is short-lived and decays approximately 
100% of the time into a b quark and a W boson, where 
W — > li^, qq'. We consider the exclusive decay mode 
t + D ^ Wb + D in which W qq' . This W decay 
mode has the largest branching ratio and it allows for 
the full reconstruction of the top quark. In this channel, 
the missing transverse energy {^t) [Si can be uniquely 
assigned to the DM particle's passage through the detec- 
tor. 

Events are collected by CDF II I!) , a general pur- 
pose detector used to study Tevatron pp collisions with 
y/s = 1.96 TeV. CDF II contains a tracking system con- 
sisting of a cylindrical open-cell drift chamber and silicon 
microstrip detectors immersed in a 1.4 T magnetic field 
parallel to the beam axis. Electromagnetic and hadronic 
calorimeters surrounding the tracking system measure 
particle energies. Drift chambers and muon scintillators 
located outside the calorimeters identify muons. We use 
a data sample corresponding to an integrated luminosity 
of 7.7±0.5 fb^i. 

We consider only those events which triggered the data 



acquisition system due to the presence of two calorime- 
ter clusters and significant ^t- We include data recorded 
between 2001 and 2010. Prior to 2007, the data acqui- 
sition system threshold was 35 GeV IJ]. After an 
upgrade to the system resulting in improved jet en- 
ergy and resolution, the requirement was lowered to 
> 30 GeV. Jets are reconstructed using the jetclu 
algorithm [2l| with a clustering radius of 0.4 in azimuth- 
pseudorapidity space ((/>, 77) [22] . Jet energies are cor- 
rected using standard techniques . Jets originating 
from b quarks are identified using a secondary-vertex- 
tagging algorithm [2^ . 

In order to retain only those events for which the trig- 
ger system is fully efficient, we select events with > 50 
GeV and three jets. Exactly one jet is identified as a 6-jet. 
We require the jet transverse energy i?;^' , to be Ei^jl > 35 
GeV, e:^ > 25 GeV, E^^^ > 15 GeV, where the jets j, 
{i — 1,2,3) are ordered by decreasing energy. We require 
that either jj or j2 have \r}\ < 0.9, and that all three jets 
have I77I < 2.4. We veto events with identified high-p^ 
electrons or muons, removing monotop events inconsis- 
tent with a hadronically-decaying top quark. 

We model the signal and background contributions to 
the selected sample using a variety of Monte Carlo (MC) 
simulation programs. In our simulation we assume a top- 
quark mass of 172.5 GeV/c^, consistent with the world's 
best determination 2^ 2g . We model monotop DM pro- 
duction in the flavor-violating process {ug — tD) with 
MADGRAPH 27|. Additional showering and hadroniza- 
tion are described by PYTHIA [2^. We have generated 11 
signal samples assuming various DM mass in steps of 5 
GeV/c^ from to 25 GeV/c^, and then in steps of 25 
GeV/c2 from 25 to 150 GeV/c^. 

The event selection described above gives a data sam- 
ple dominated by QCD multijet events, where the false 
arises from the mismeasurement of jet energy. Sim- 
ulation of this background is prohibitive due to the high 
production rate and large theoretical uncertainties. In- 
stead, we use a method which relies on data and is 
based on a recently improved Tag Rate Matrix (TRM) 
method [2^ . The TRM method utilizes an estimate of 
the probability for QCD multijet events to have tagged 
jets. The probability is derived in a control region dom- 
inated by QCD multijet events. This probability is ap- 
plied as a per-event weight to all events meeting our anal- 
ysis selections excluding the 5-jet requirement. From this 
sample of weighted events, we subtract the expected elec- 
troweak components (as modeled by applying the same 
TRM probability to simulated samples). The resulting 
events form our model of the QCD multijet component 
of the analysis data sample. 

We model other physics with samples generated by 
MC programs. Diboson and tt production are gener- 
ated by PYTHIA and normalized to the next-to-leading 
order (NLO) cross section predicted using the mcfm 
program [sij and the approximate next-to-next-to- 
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leading order cross section [32'| , respectively. The pro- 
duction of W/Z plus light flavor and heavy flavor (HF) 
jets are simulated by ALPGEN 33| with showering and 
hadronization performed by PYTHIA and normalized to 
NLO cross sections. Single top, both s— and t- channel 
production, are modeled using MADGRAPH with PYTHiA 
and normalized to NLO cross sections 3J, |35 [ . 

The light flavor jets misidentified as 6-jets by 
the secondary-vertex-tagging algorithm are labeled as 
mistags. A data-driven method is used to estimate the 
mistag rate for the tagging algorithm [2^ . We apply the 
mistag rate to the MC events with light flavor jets to 
estimate the mistag contribution. 

Figure[l] shows the distribution in a control region 
for events which pass our signal selection but have an 
identified high-p^ electron or muon. 
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FIG. 1: The distribution in a control region requiring an 
identified high-pT lepton; the gray area represents the statis- 
tical and systematic uncertainties on the background model. 
The lower panel displays the difference between the data and 
the expected backgrounds (Exp) divided by the latter. 

After the selection described above, we are left with 
6471 data events. We expect that approximately 70% 
of these events come from QCD multijet production. In 
order to further suppress the QCD contamination and 
the other SM backgrounds, we require the azimuthal dis- 
tances between the and j2, A0(^t, j2) > 0.7, as the 
in QCD multijet background tends to align to the jet 
with less measured energy. We also require the invariant 
mass of the three jets to be consistent with the recon- 
structed top-quark mass, 110 < rrijjj < 200 GeV/c^, 
large significance i^T/\/Yl > S.SVCeV, where 
J2 Et is the scalar sum of transverse energy deposited 
in the calorimeter) and Eji > 25 GeV. All selections 
have been chosen to optimize the significance S/y^S + B, 
where S and B are the expected number of signal and 



backgrounds events, respectively. Table|T]shows the num- 
ber of events in the signal region for the data, the number 
of events for SM backgrounds, and the expected signal 
assuming different values of the DM particle's mass. The 
events that fail these signal-region selections are used to 
form a control region that is used to validate the back- 
ground models, as well as to determine the normalization 
of the QCD multijet background. 



TABLE I; Number of expected signal and background events 
compared to data in the signal region. The expected signals, 
assuming different values for the mass of the DM particle, are 
also presented. The errors include statistical and systematic 
uncertainties. 



Processes 


Events 


pp^t + D 




mo ^20 GeV/c^ 


2116.9 ± 121.4 


mc = 75 GeV/c^ 


232.3 ± 22.9 


mo = 100 GeV/c^ 


129.8 ± 12.5 


mo = 125 GeV/c^ 


94.5 ± 9.3 


tt 


182.8 ± 20.2 


Single top 


24.3 ± 4.5 


Diboson 


15.7 ± 2.7 


W/Z+RF 


130.5 ± 33.8 


Mistag 


96.9 ± 39.4 


QCD multijet 


210.2 ± 54.5 


Total background 


660.2 ± 78.1 


Data 


592 



We consider several systematic uncertainties affecting 
the sensitivity of this search. The dominant systematic 
sources are the uncertainties on multijet normalization 
(25.5%), the mistag rate (16.6%) and the background 
cross sections (6.5% — 30%). We also consider uncertain- 
ties from the jet energy scale [ll] (2.8% - 10.7%), the 
luminosity measurement [111 (6%), parton density func- 
tions (2%), lepton veto (2%), 6-tagging efhciency (5.2%), 
trigger efficiency (0.4% — 0.9%), and from the initial-state 
and final-state radiation (4%). We also assign systematic 
uncertainties, based on the variation in the shape of the 
distribution of kinematic quantities, under a itlcr varia- 
tion of the jet energy scale and the uncertainty on the 
efficiency of the data acquisition system. 

The is chosen to discriminate the signals from the 
backgrounds. The distribution due to a DM parti- 
cle of mass of 125 GeV/c^ and the SM backgrounds are 
shown in Fig. [51 The signal is expected to contribute 
significantly at high values of ^p. We find no significant 
excess of signal- like events in the data analyzed, and thus 
proceed to set 95% confidence level (C.L.) upper limits 
on the monotop DM production cross section. The limits 
are calculated with the distribution as the shape dis- 
criminant using a Bayesian maximum likelihood method 
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assuming a flat prior for the signal cross section 37[ . We 
treat systematic uncertainties using a Bayesian marginal 
likelihood method. Figure[3] shows the calculated upper 
limits on the monotop cross section as a function of the 
mass of the DM candidate compared to the theoretical 
predictions. 
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FIG. 2: The distribution in the signal region. The data 
is compared to the sum of the SM contributions. The distri- 
bution of signal events with a DM mass of 125 GeV/c^ is also 
shown. 
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FIG. 3: Exclusion curve of the monotop cross section as a 
function of the mass of DM particle 

In conclusion, we have performed the first search for 
the production of DM in association with a single top 
quark at hadron colliders. In an analysis of 7.7 fb~^ of 
CDF II data we have found that the observed data is 
consistent with the expectation from SM backgrounds. 
We set 95% CL. upper limits on the cross section of 
pp D + t as a. function of the DM mass in the range 
of — 150 GeV/c^. Future searches for new physics in 
monotop final states can probe resonant production of 



top quarks and DM candidates with exotic mediators. 
While these processes are predicted to have low produc- 
tion rates (making them difficult to probe with Tevatron 
data), they are expected to be within the reach of LHC 
experiments with sufficient data. 
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